Atomic carbon chains as spin-transmitters: an A b initio transport study 
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An atomic carbon chain joining two graphene flakes was recently realized in a ground-breaking 
experiment by Jin et at, Phys. Rev. Lett. 102, 205501 (2009). We present ab initio results 
for the electron transport properties of such chains and demonstrate complete spin-polarization 
of the transmission in large energy ranges. The effect is due to the spin-polarized zig-zag edge 
terminating each graphene flake causing a spin-splitting of the graphene n z bands, and the chain 
states. Transmission occurs when the graphene 7r-states resonate with similar states in the strongly 
hybridized edges and chain. This effect should in general hold for any ^-conjugated molecules 
bridging the zig-zag edges of graphene electrodes. The polarization of the transmission can be 
controlled by chemically or mechanically modifying the molecule, or by applying an electrical gate. 
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Introduction. Carbon based materials have assumed a 
central role in the development of nanoelectronics 
The spectacular electronic properties of Carbon nan- 
otubes and 2D graphene sheets are exploited in many 
device proposals possibly paving the way to an era of 
carbon-based electronics. Further, the magnetic proper- 
ties of these structures make them promising candidates 
in spintronics and numerous suggestions have been made 
for carbon based spin-filtering systems [2Hl4j. 

The fabrication of a field effect transistor by carv- 
ing a constriction in graphene has triggered a mas- 
sive interest in such geometries. Graphene nanoribbons 
forming the constrictions have been studied intensively. 
Already some time ago Nakada et al. [§] showed that 
these ribbons have peculiar edge states which depend on 
the edge geometry as well as the ribbon width. Zig-zag 
edges are particularly intriguing as they support local- 
ized spin-polarized states. For example, Muhoz-Rojas et 
al. predict a giant magnetoresistance in a device based 
on zig-zag ribbons connected to metal electrodes [l2j ]. 
The thinnest possible constriction, the monoatomic lin- 
ear chain, was already several years ago proposed theo- 
retically as an ideal component in molecular devices 
llij ]. Until recently these chains have been studied only 
in very few experiments due to the daunting challenges 
in their fabrication 1^- 21 j . 

Recently Jin et al 



22j realized quite stable freestand- 



ing carbon chains connected to two graphene flakes. By 
employing energetic electron irradiation inside a trans- 
mission electron microscope two holes were created in a 
large graphene sample. The resulting ribbon separating 
the two holes was carefully thinned by lower energy elec- 
trons until only the monoatomic chain remained. This 
structure is thus both the realization of a pure carbon 
constriction of the thinnest possible kind, as well as an 
atomic chain contacted by carbon leads instead of the 
usual case of metallic contacts. Chuvilin et al. have re- 
ported similar results [3] ■ Jin et al. also performed den- 
sity functional theory (DFT) studies showing that the 



zig-zag edge is the energetically most favorable structure 
for connecting the chain, in agreement with the experi- 
mental observations. Recently Chen et al. 24| reported 
calculations of the electronic transmission through this 
system based on a tight-binding 7r-orbital model; the 
dominant features were asymmetric resonant peaks close 
to the Fermi level. 




FIG. 1: (Color online) The C-9 system with 9 carbon atoms 
in the chain. The shaded areas marked L and R are the 
electrodes, and C is the central region. Periodic boundary 
conditions are imposed in the transverse, x, direction. 



In this Letter, we describe the spin-dependent elec- 
tronic transport properties of carbon chains connect- 
ing two graphene sheets at the zig-zag edges. From ab 
initio calculations we find a strong spin-polarization of 
the transmission caused by the intrinsic properties of the 
zig-zag edge: it spin-polarizes and introduces a spin split- 
ting of the graphene it bands as well as the hybridized 
edge-chain states. The peaks in the transmission aris- 
ing as graphene states are in resonance with the strongly 
hybridized edge and chain states are thus separated in 
energy resulting in the polarization of the current. We 
argue that this is a generic feature of this class of systems 
making them strong candidates for future spintronic de- 
vices. 

Method and Systems— We use ab initio pseudopotential 
spin-polarized DFT as implemented in the Siesta code 
25j to obtain the electronic structure and relaxed atomic 
positions. Our spin transport calculations are based on 
the nonequilibrium Green's function method [26] as im- 
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FIG. 2: (Color online) (Top): The electron transmission of 
the C-8 and C-9 systems. (Middle): PDOS of chain atom n z - 
orbitals. (Bottom): The spin-polarization of the transmission 
shown above given in terms of TSP. 



plemented in the Transiesta[27| code, extended to spin- 
polarized systems [40|. 

The model system is inspired by that of Ref. 22j and 
is shown in Fig. [TJ in the case of N = 9 carbon atoms 
in the chain, denoted C-9. The system is divided into 
left (L) and right (R) electrodes, and a central region C 
(27j . The electrodes are semi-infinite and we employ peri- 
odic boundary conditions in the transverse, x, direction. 
Most importantly, the two graphene sheets are termi- 
nated by zig-zag edges. Initially, we consider them to be 
H-passivated. We use a fixed C-C distance of 1.3 A in 
the chain. We find no qualitative differences in the trans- 
mission from using fully relaxed odd N structures 41 1. 



For even N a careful relaxation migh result in a dimeri- 
sation in the chain as found by Khoo et al. 28|. To test 



whether this affects the transmission, we have computed 
the transmission with fixed single/triple bond lengths of 
1.26/1.36 A in the chain of the C-8 system, as discussed 
below. 

Results. The spin-resolved transmission is shown in 
Fig. Eftop), for two systems with different chain lengths. 
Both cases show large peaks away from the Fermi level 
Ep, for both spin channels and smaller peaks closer to 
Ep. All peaks have asymmetric line-shapes in agreement 
with Chen et al. [24j. At Ep the transmission is zero. 
Outside the energy window |E| < 1 eV shown in Fig. 
[21 large peaks are found repeatedly at higher and lower 
energies, but with much less pronounced spin-splitting. 
We focus here on the range |E| < 1 eV around the in- 
trinsic Fermi-level, which can be modified by chemical 
doping [29j], or by a gate voltage [l5j]. In this range we 
obtain very similar transmissions for C-7 and C-9, and 
for C-8 and C-10, respectively. Thus, a clear odd/even 



effect is present as in other carbon chain systems 16|, [17 1 . 
The introduction of a rigid dimerisation for the C-8 sys- 
tem causes minor relative shifts in the resonance peaks 
but the qualitative picture is unaltered (these results are 
available in the Supplementary material SI). We stress 
that a detailed study of odd/even effects is not our pri- 
mary concern here, and that the comparisons discussed 
here only serve to demonstrate that the spin-polarization 
persists in both cases. 

The total transmission is spin-polarized as the peaks 
of each spin channel are shifted relative to each other. 
This is illustrated in Fig. [2] (bottom), where we plot the 
spin-polarization of the transmission (TSP) given by 



TSP = 



Min 



Maj 



(1) 



Significant energy windows above and below Ep dis- 
play complete polarization, |TSP| ~ 1. The polarization 
changes sign around Ep being predominantly minority 
(majority) spin below (above) Ep. The majority (minor- 
ity) spin is defined as having the highest (lowest) electron 
occupation. 

We next examine the features in the transmission in 
more detail. An infinite monatomic carbon chain has 
two spin-degenerate transmitting channels; one channel 
consists of the 7r z -orbitals and the other one of the tt x - 
orbitals not participating in the sp hybridization in the 
chain. Since transport takes place via the 7r z -orbitals 
in the graphene sheets and since they are orthogonal to 
the 7r x -orbitals, transmission will occur only through the 
7r z -orbitals of the system. This explains why the peak 
height is close to unity for each spin channel, as seen in 
the transmission spectra of the total system in Fig. [2] 
(top). 

Transmission occurs when graphene states are in res- 
onance with the states of the atomic chain. This is seen 
from the projected density of states (PDOS) of chain 
atom 7r z -orbitals plotted in Fig. [5] (middle), which shows 
a clear correlation with the transmission. The height 
of the peaks in transmission reflects the linear DOS of 
graphene. The chain states differ from those of a corre- 
sponding free chain as the chains in our case are found 
to strongly hybridize with the graphene edges. This is 
also reflected in the magnetic moment of the contact 
edge atom being lower by roughly 0.2 /is compared to 
the edge atom furthest away from the chain. The spin- 
splitting of the transmission curves originates from the 
spin-polarization of the two zig-zag edges. The energetics 
favors a ferromagnetic chain-mediated coupling between 
the edges, e.g., E G {AFM) - E G (FM) = 0.29 eV for a 
nine-atom chain where Ec(AFM) is the total energy of 
the AFM solution (See also SI). A (weak) external mag- 
netic field is needed to define a preferred spin orientation 
in an experiment. In our calculations we consider only 
the ferromagnetic case. The magnetic moment of the 
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FIG. 3: (Color online) The real part of the minority spin 
eigenstates of the C-9 system (a,c) and transmission eigen- 
channels (b,d) at the resonances indicated by 1 and 2 in Fig. 
[21 respectively. 



chain, which is N dependent, almost entirely occurs in 
the 7Ta;-orbitals and has little influence on the transport. 

The resonating states are investigated further in Fig. 
[3l using the minority spin states of the C-9 system as 
an example. In (a) and (c) are shown the real part of 
the eigenstates of the total system with energies corre- 
sponding to the peaks indicated in Fig. [2] by 1 and 2, 
respectively. The eigenstates show a strong hybridiza- 
tion between the edges and the chain. The transmission 
eigenchannels at the resonances are shown in (b) and 
(d), respectively. Except for a phase shift, the matching 
with the eigenstates is clear. At energies corresponding 
to zero transmission we find no eigenstates involving the 
7T z -orbitals. Calculations on systems with wider leads 
confirm that hybridization only involves the part of the 
edge close to the chain. The same analysis can be carried 
out for the other resonances, and for both spins. 

We have tested the robustness of the spin-polarization 
effect in various ways (details available in SI): (i) apply- 
ing a bias voltage of 0.5 eV (ii) gating with a lateral (z) 
field of 0.25 V/A ; (iii) adding a charge of -0.5 e; and (iv) 
doubling the electrode width which reduces the relative 
weight of the chains along the edge. Also, the electrodes 
in a real sample are not semi-infinite sheets, but rather 
large ribbons or flakes. We therefore also examined 9- 




FIG. 4: (top) The transmission for a zig-zag shaped carbon 
chain which is partly H terminated, (bottom) The transmis- 
sion of the C-9 system stretched by 1 A. The results in the 
unstretched case are shown for comparison (dotted line). 
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FIG. 5: TSP for a 9 atom chain connected to partly H- 
passivated(top) and unpassivated (bottom) edges. 



atom wide fully H-passivated semi-conducting armchair 
ribbon electrodes. In that case the details in the trans- 
mission close to Ep are influenced by the band gaps in 
the leads, but the spin-polarization of the resonances is 
still present. All these tests confirm the qualitative pic- 
ture described above. 

The fact that the spin-polarization of the transmission 
is caused by an intrinsic property of the zig-zag edge sug- 
gests that the high TSP values might be a generic prop- 
erty of these systems. The spin-polarization effect should 
then occur also for all carbon molecules which make an 
sp 2 bond to the zig-zag edge, and conduct via the ir z - 
orbitals. In order to test this hypothesis we have per- 
formed calculations for a system where the linear chain is 
replaced with a partly H-passivated zig-zag shaped chain 
as shown in the inset in Fig. |4j upper panel. Again, the 
splitting of the large resonances appears in the same man- 
ner as for the linear chain, however even more strongly 
in this case. Yet another possibility for tuning the TSP 
is to mechanically change the electrode distances which 
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induces changes in the chain states. This is shown in Fig. 
IH lower panel, where the C-9 system has been stretched 
by 1 A. 

Finally, we comment on the influence of defects and 
edge passivation. Edge disorder has been shown to af- 
fect the transport properties of graphene ribbons 
However, the crucial point in our case is whether the mag- 
netisation is quenched or not. The chain acts as a defect 
which is reflected in the mentioned significant lowering of 
the magnetic moment of the contact atom. However, we 
find that the magnetisation along the - relatively short 
- edge is rapidly restored away from the contact atom, 
in agreement with other studies 35j, |36[ . We thus infer 



that the effect is expected even in the presence of nearby 
defects. The extent of edge passivation in actual samples 
is uncertain 37 , 3^1 ■ In the results presented so far, we 



have considered full H passivation. In Fig. [S]we show the 
TSP for the C-9 system in the case where two H atoms 
has been removed from one edge so half of the edge atoms 
are passivated, as well as the case of no passivation at all. 
The TSP is seen to be largely intact below Ep, but re- 
duced above Ep compared to the case of full passivation 
shown in Fig. [5] The TSP above Ep can, however, be 
tuned as explained above. 

Conclusion. We have presented ab initio results for 
the electron transport properties of atomic carbon chains 
between zig-zag edges of graphene electrodes. We find 
a complete spin-polarization of the electron transmis- 
sion in large energy ranges around the Fermi level owing 
to the intrinsic spin-polarization of the zig-zag edges of 
graphene. Transmission occurs as graphene states are in 
resonance with the strongly hybridized edge and chain 
7r 2 -orbital states. The spin-polarization of the trans- 
mission can be controlled by manipulating the chemical 
structure of the chain or a similar carbon molecule, ap- 
plying mechanical strain, or changing the position of the 
Fermi-level by doping or gating. 

Acknowledgments. Computational resources were pro- 
vided by the Danish Center for Scientific Computations 
(DCSC). APJ is grateful to the FiDiPro program of 
Academy of Finland. 



* Corresponding author: joachim.fuerst@nanotech.dtu.dk 

[1] P. Avouris, Z. Chen, and V. Perebeinos, Nature Nan- 
otech. 2, 605 (2007). 

[2] J. Ferrer and V. M. Garcia-Suarez, Journal Of Materials 
Chemistry 19, 1696 (2009). 

[3] N. Tombros, C. Jozsa, M. Popinciuc, H. T. Jonkman, 
and B. J. van Wees, Nature 448, 571 (2007). 

[4] J. Fiirst, J. Pedersen, C. Flindt, N. A. Mortensen, 
M. Brandbyge, T. Pedersen, and A. P. Jauho, New Jour- 
nal of Physics 11, 115117 (2009). 

[5] M. Vanevic, V. M. Stojanovic, and M. Kindermann, 
Phys. Rev. B 80, 045410 (2009). 



[6 
[7 
[8 
[9 
[10 

[11 

[12 

[13 

[14; 

[15 
[16 

[17; 

[18 
[19 

[20 
[21 
[22 
[23 
[24 
[25 

[26 

[27 
[28 
[29 
[30 

[31 
[32 
[33 
[34 
[35 
[36 



D. Bahamon, A. Pereira, and P. Schulz, Phys. Rev. B 89, 
125414 (2009). 

O. V. Yazyev and A. Pasquarello, Phys. Rev. B 80, 
035408 (2009). 

M. Fujita, K. Wakabayashi, K. Nakada, and K. Kusak- 

abe, J. Phys. Soc. Jap. 65, 1920 (1996). 

K. Nakada, M. Fujita, G. Dresselhaus, and M. S. Dres- 

selhaus, Phys. Rev. B 54, 17954 (1996). 

W. L. Wang, S. Meng, and E. Kaxiras, Nano Letters 8, 

241 (2008). 

O. V. Yazyev, Nano Letters 8, 1011 (2008). 

F. Munoz-Rojas, J. Fernandez-Rossier, and J. J. Pala- 

cios, Phys. Rev. Lett. 102, 136810 (2009). 

J. A. Fiirst, M. Brandbyge, A. P. Jauho, and K. Stokbro, 

Phys. Rev. B 78, 195405 (2008). 

M. Koleini, M. Paulsson, and M. Brandbyge, Phys. Rev. 
Lett. 98, 197202 (2007). 

A. K. Geim and K. S. Novoselov, Nature Mat. 6, 183 
(2007). 

N. D. Lang and P. Avouris, Phys. Rev. Lett. 81, 3515 
(1998). 

B. Larade, J. Taylor, H. Mehrez, and H. Guo, Phys. Rev. 
B 64, 075420 (2001). 

S. Tongay, R. T. Senger, S. Dag, and S. Ciraci, Phys. 
Rev. Lett. 93, 136404 (2004). 

R. J. Lagow, J. J. Kampa, H. C. Wei, S. L. Battle, J. W. 
Genge, D. A. Laude, C. J. Harper, R. Bau, R. C. Stevens, 
J. F. Haw, et al., Science 267, 362 (1995). 
X. L. Zhao, Y. Ando, Y. Liu, M. Jinno, and T. Suzuki, 
Phys. Rev. Lett. 90, 187401 (2003). 

T. D. Yuzvinsky, W. Mickelson, S. Aloni, G. E. Begtrup, 
A. Kis, and A. Zettl, Nano Letters 6, 2718 (2006). 

C. H. Jin, H. P. Lan, L. M. Peng, K. Suenaga, and 
S. Iijima, Phys. Rev. Lett. 102 (2009). 

A. Chuvilin, J. C. Meyer, G. Algara-Siller, and U. Kaiser, 
New Journal of Physics 11, 083019 (2009). 
W. Chen, A. Andreev, and G. Bertsch, Phys. Rev. B 80, 
085410 (2009). 

J. M. Soler, E. Artacho, J. D. Gale, A. Garcia, J. Jun- 
quera, P. Ordejon, and D. Sanchez-Portal, Journal Of 
Physics-Condensed Matter 14, 2745 (2002). 
H. Haug and A. P. Jauho, Quantum kinetics and optics 
of Semiconductors, 2nd revised edition (Springer Verlag, 
Berlin, 2008). 

M. Brandbyge, J.-L. Mozos, P. Ordejon, J. Taylor, and 

K. Stokbro, Phys. Rev. B 65, 165401 (2002). 

K. H. Khoo, J. B. Neaton, Y. W. Son, M. L. Cohen, and 

S. G. Louie, Nano Letters 8, 2900 (2008). 

K. Novoselov and A. Geim, Materials Technology 22, 178 

(2007). 

M. Y. Han, B. Ozyilmaz, Y. B. Zhang, and P. Kim, Phys. 
Rev. Lett. 98, 206805 (2007). 

D. Gunlycke, D. A. Areshkin, and C. T. White, Appl. 
Phys. Lett. 90, 142104 (2007). 

D. A. Areshkin, D. Gunlycke, and C. T. White, Nano 
Letters 7, 204 (2007). 

M. Evaldsson, I. V. Zozoulenko, H. Y. Xu, and 
T. Heinzel, Phys. Rev. B 78, 161407 (2008). 
H. Y. Xu, T. Heinzel, and I. V. Zozoulenko, Phys. Rev. 
B 80, 045308 (2009). 

D. E. Jiang, B. G. Sumpter, and S. Dai, J. Chem. Phys. 
127, 124703 (2007). 

M. Wimmer, I. Adagideli, S. Berber, D. Tomanek, and 
K. Richter, Phys. Rev. Lett. 100, 177207 (2008). 



[37] C. O. Girit, J. C. Meyer, R. Erni, M. D. Rossell, 
C. Kisielowski, L. Yang, C.-H. Park, M. F. Crommie, 
M. L. Cohen, S. G. Louie, et al., Science 323, 1705 
(2009). 

[38] T. Wassmann, A. P. Seitsonen, A. M. Saitta, M. Lazzeri, 
and F. Mauri, Phys. Rev. Lett. 101, 096402 (2008). 

[39] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

[40] We employ GGA PBE for exchange-correlation [3{|, 
mesh cutoff of 200 Ry, and Monkhorst-Pack k-sampling 
of (4,1,10). The transmission is averaged over 100 Ap- 
points. The basis set size is single-^ in transport calcula- 
tions, and we checked that larger double-^ does not have 
qualitative influence on the transmission. 

[41] We have performed relaxation with a double-^ basis set 
on the C-9 system using the conjugant gradient method 
and a force tolerance of 0.01 eV/A. Only the chain atoms 
and the two closest rows of atoms on each side were al- 
lowed to move. We find bonding lengths in agreement 
within 0.05 A of Ref. 



6 

Supplementary material 



Lateral gate 

c-9 




Charge -0.5e 

C-9 




E-E F (eV) 



Bias of 0.5 V 

C-9 



Majority spin 

Minority spin 




Double electrode width 

C-9 




-1 -0.5 0.5 1 



E-E F (eV) 



FIG. 6: The electron transmission for the C-9 system when (top left) applying a bias voltage of 0.5 eV, (top right) gating 
with a lateral (z) field of 0.25 V/A , (bottom left) adding a charge of -0.5e, (bottom right) doubling the electrode width which 
reduces the relative weight of the chains along the edge. 



(a) 



(b) 



FIG. 7: The difference in majority and minority spin densities for the C-9 system. Bronze (turquoise) indicates surplus 
of majority (minority) spin, (a): The FM solution between edges. The spin resides on the 7r z - orbitals at the edges and 
predominantly on the n x - orbitals in the chain, (b): An AFM solution obtained by initially breaking the FM symmetry by 
hand. Notice, that the chain is largely spin unpolarized. We obtain E^~ 9 {AFM) - E^~ 9 (FM) = 0.29 eV and E%~ a {AFM) - 
Eq~ 8 (FAI) = 0.31 eV. The connection of the two edges via the chain thus seem to favor a FM solution. 
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FIG. 8: The transmission for a 9 atom chain connected to fully H-passivated armchair semiconducting ribbon electrodes with 
a width of 9 rows of atoms. 
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FIG. 9: The transmission (top) and TSP (bottom) for the C-8 system where a dimerisation has been inserted by hand with 
single/triple-bond lengths of 1.26/1.36 A. 



